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Abstract 

The purpose of this paper is to provide an introduction to the physics 
of scattering theory, to define the dosimetric concept of linear energy 
f— | ■ transfer in terms of scattering theory, and to provide an introduction to 

D ' the concepts underlying Monte Carlo simulations. 

O ■ 1 Introduction 

£>y Whilst the absorbed dose deposited by energetic particles ('radiation') in con- 

4h ' densed matter is simply the energy per unit mass delivered to the medium, the 

equivalent dose is a function of both the energy per unit mass, and the spa- 
tial distribution with which that energy is deposited. The greater the spatial 
concentration of deposited energy, the greater the equivalent dose. The average 
rate at which a type of radiation locally deposits energy in a medium, per unit 
distance of track length, is called the linear energy transfer (LET) of that radi- 
ation in that medium. Radiation with higher linear energy transfer deposits a 
greater equivalent dose in a medium. 

The LET is closely related to the stopping power of the medium, which 
■ is the energy lost by an incoming particle, whether it is deposited locally, or 

OO ! transformed into other radiation, such as Bremsstrahlung photons, or the en- 

ergy of secondary electrons, (so-called 'delta-rays'). Linear energy transfer and 
stopping power generally coincide in the case of heavy charged particles, but 
the linear energy transfer of beta-radiation (i.e, incoming electrons) does not 
K% ' include the energy transformed into Bremsstrahlung photons or delta-rays. 

^ ■ The first purpose of this paper is to define linear energy transfer in terms 

of the concepts used in scattering theory, and, in particular, the concepts of 
quantum scattering theory. An exposition will then be provided of the principles 
underlying Monte Carlo simulations. 



2 Scattering 

Scattering theory is the branch of physics which describes the collision processes 
between particles. These are typically considered to be processes in which an 
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incoming particle, or beam of such particles, interacts with a target particle, or 
collection of particles. After an individual interaction, the energy and momen- 
tum of both the incoming particle and the target particle can be altered. A 
change of momentum includes both a change in the direction of travel, and, in 
the case of particle with non-zero mass, a change in the speed of travel. The 
target particle is often considered to be a composite particle, and as such, both 
its energy of motion, and its internal energy state, can change as a result of the 
collision. 

Consider first the scattering of a particle beam in terms of classical physics. 
In the case of a beam of particles with non-zero mass, the beam has a flux F 
defined as 



F = n i v l , 

where rii is the number of particles per unit volume in the beam, and vi is 
the average velocity of the particles in the beam. The flux has dimensions 
of Area" 1 Time" 1 . It can be equivalently defined as the number of particles 
passing though a unit cross-sectional area of the beam, per unit time. The flux 
should not be confused with the fluence $ of a particle beam, which is the 
number of particles passing through a unit cross-sectional area of the beam, 
over the time duration considered. Hence, the fluence is the definite integral of 
the flux over time, $ = f t F(t) dr, and the flux is the time derivative of the 

fluence, F = d<fr/dt. The fluence has dimensions of Area" 1 . 

A particle detector placed at an angle = (9, <f) with respect to the direction 
of the incoming beam, will detect particles at the rate (Bohm 1979, p312): 

N/At = a(Q)N T n % Vi 
= a(Q)N T F . 

Nt is the number of particles in the target, and the constant of proportionality 
cr(f2) is called the differential cross-sectionj in the direction fi. The differential 
cross-section can clearly be expressed as 

o-(fi) - 1 N/AT 



or, in words, 



a(Q) 



Nt riiVi 



1 Number of particles scattered into Q per unit time 
Nt Incident Flux 



The dimensions of the differential cross-section are Area ■ Solid Angle" 1 . 

Integrating the differential cross-section over all possible directions gives the 
total scattering cross-section a: 



1 In the physics literature, the differential cross-section is usually denoted as da, or da/dQ.. 
The expression for da typically contains a dependence upon both the direction and magnitude 
of the outgoing momentum, whilst the expression for da/dQ contains only a dependence upon 
direction, and is obtained by integrating over all the outgoing momenta in that direction. We 
shall refrain from such mathematically questionable notation in this paper. 
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(j = J a(fl)dn , 

where dfl — sin 9d(j) e?#0 The dimensions of the total cross-section are simply 
Area. One can think of the cross-section as effectively the area which the target 
presents to the incoming beam. 

In terms of the quantum mechanics of a scattering process, an incoming 
particle is represented by a quantum state, or 'wave-function' ^j, the interaction 
process is represented by a scattering operator S, the final outgoing particle state 
is represented by another wave- function ^ f, and the objective is to calculate the 
transition probability |(^E'/|iS'\E'i}| 2 between an incoming state and an outgoing 
statejfl In many physically relevant situations, the incoming state has a specific 
energy Ei and momentum ki , and each possible outgoing state also has a specific 
energy Ef and momentum kf. An outgoing state with a specific momentum kf, 
also has a specific direction Q associated with it. 

The square-amplitude \^i\ 2 of the incoming particle's wave-function gives 
the position probability density for the particle. In quantum mechanical terms 
the incoming flux is then the product of the probability density per unit 
volume, with the velocity of the particle. Equivalcntly, the quantum mechanical 
flux is the probability density per unit area, per unit time. Hence, in quantum 
mechanical terms even an individual particle has a flux associated with it. Given 
the transition probability | / {S^i) | 2 , there is also an associated transition rate, 
or transition probability per unit time. In the simplest quantum mechanical 
terms, cross-sections can be defined as follows (Bohm 1979, p314): 

Transition probability 

a = , 

Incident probability per unit area 

or, equivalently, 

Transition probability, per unit time 

a = ; — ; ; ; , 

Incident probability per unit area, per unit time 
which itself is equivalent to 

Transition rate 
Incident Flux 

The quantum mechanical cross-sections used in practice provide an econom- 
ical way of bundling the transition probabilities between classes of quantum 
states. For example, the differential cross-section <r(E, il) is proportional to 
the probability of a transition from any incoming state of energy E to any 
outgoing state ^ f in which the momentum vector k f points in the direction of 
fi. As before, integrating a differential cross-section over all possible directions 
then gives a total scattering cross-section <j(E): 

2 In terms of differential geometry, dQ is the volume form on the 2-sphere, expressible as 
dQ = sin6W0 A dd, where A is the antisymmetric tensor product. 
3 (-|-) here is the inner product on the space of quantum states. 
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a(E) = / cr(E,n)dn 



This total cross-section is proportional to the probability of an incoming state 
of energy E interacting with the target particle. 

Now, an elastic collision is defined to be a collision in which the outgoing 
particles are the same as the incoming particles, and the total kinetic energy of 
the particles before the collision is the same as the total kinetic energy of the 
particles after the collision. Interactions in which this does not occur are called 
inelastic collisions. For example, an inelastic collision occurs when an incoming 
particle interacts with an atom, and some of the energy from the incoming 
particle excites or ionises the atom. Alternatively, the incoming particle can 
be absorbed by the target, raising the target particle into an excited state, 
which may then subsequently release various decay products. An example of 
this is the absorption of an incoming neutron by an atomic nucleus. Hence, 
there are elastic scattering cross-sections, inelastic scattering cross-sections, and 
absorption cross-sections, all of which are functions of the energy of the incoming 
particle(s). 

Given the inelastic differential cross-section <T^(f2) for an interaction in which 
the incoming projectile transfers AE^ of energy to a target atom, and the 
projectile scatters with a momentum kf in direction f2, the total inelastic cross- 
section for all interactions in which the projectile transfers AE^ of energy to a 
target atom is: 



Given N target atoms per unit volume in a condensed matter medium, the 
product Na^ is the probability per unit path-length of the incoming particle 
track, that the projectile will lose AE^ of energy. Hence, the local energy 
transfer (linear energy transfer) of the particle is given by the expression: 



4 Given the probability pi for each possible value ct^ of a quantity A, the expectation value 
(A) of the quantity is defined to be the probability-weighted sum of those values, (A) = 
J2i ViO-i- 





To be more precise, this is the quantum mechanical expectation valuc@ for the 
energy loss per unit length: 
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A special case of this equation is Bethe's formula for the local energy loss, per 
unit track-length, of a swiftly moving charged particle heavier than an electron. 
Bethe's formula applies, for example, to energetic protons, alpha particles, and 
heavy ions. For swift charged particles heavier than an electron, the inelastic 
cross-sections are sharply peaked towards the forward direction. If we take the 
forward direction to correspond to a polar angle of 9 = 0, such particles scatter 
into a narrow cone-shaped region circumscribed by a small polar angle O, and 
all azimuthal angles <\> € [0, 2tt). Bethe's formula then arises from assuming that 
the cross-section 



Assuming the projectile has a non-relativistic speed, Bethe's formula for linear 
energy transfer is: 



where M is the mass of projectile, v is its speed, ze is its charge, m e is the 
electron mass, Z is the proton number of the target atoms, N is the number of 
target atoms per unit volume, and / is a parameter, called the 'mean excitation 
potential', characterising only the target atoms. Because the mass M of the 
projectile only appears in the natural logarithm, the charge of the incoming 
projectile has a greater effect upon linear energy transfer than the mass of the 
incoming projectile. 

3 Monte Carlo Simulations 

Monte Carlo simulations are computer simulations which use random numbers. 
Given the inherently probabilistic nature of particle interactions, Monte Carlo 
simulations can be used to simulate particle tracks, and their local energy de- 
position in a chosen medium. Monte Carlo simulations can thereby be used to 
calculate absorbed dose and equivalent dose. The cross-sections obtained from 
quantum mechanical and quantum field theoretical scattering theory provide 
the probability distributions used in such Monte Carlo simulations. 

A computer can generate random numbers with any specified probability 
distribution. In other words, a computer can generate numbers with the rel- 
ative frequencies specified by any probability distribution. The ability of a 
computer to generate such random numbers is dependent upon the ability of 
a computer to act as a 'pseudo-random' number generator. The computer im- 
plements an algorithm which generates random numbers in the interval (0, 1) 




can be approximated by 
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with an apparently uniform distribution. The algorithm generates each number 
from its predecessor according to a deterministic rule, hence the randomness is 
purely one of appearance. The starting point of the pseudo-random generator 
is called the 'seed', and if one starts successive runs of the program with the 
same seed, then exactly the same sequence of numbers will be generated, in the 
same order, in each run. If one wishes to avoid generating the same sequence of 
numbers on successive runs, then one simply changes the seed of the generator. 
The important point, however, is that the numbers are, to a specified level of 
approximation, uniformly distributed in the interval (0, 1). 

Given the generation of random numbers with a uniform distribution in 
(0, 1), random numbers can be generated in any other interval, with any specified 
distribution. The method of doing so exploits the fact that because the definite 
integral of any probability distribution p(x) must equal one, 



the indefinite integral F(x), which uses the upper limit of integration as the 
dependent variable, 



must assume values in the interval between and 1, beginning at the value of 
when x = a, and increasing until it reaches the value of 1 when x = b. 

Given a probability distribution p(x) over a domain (a, b) C R, the selection 
of a random number p, between and 1, selects a value ( G (a, b) such that 



When many values in (a, b) are selected by this means, they will approximate 
the distribution specified by p(x). The selection of random values according to 
a specified probability distribution, in this manner, is referred to as 'sampling 
from a random distribution.' 

In the Monte Carlo simulation of a particle track through a medium, one 
typically generates random numbers for (i) the distance travelled by the particle 
before its next collision; (ii) whether the next collision is elastic or inelastic; (iii) 
the type of inelastic collision in such an event; (iv) the energy transferred in the 
event of an inelastic collision; and (v) the direction with which the projectile 
scatters as a consequence of the collision. Let us look at each of these in turn. 

To calculate the distance travelled I by a particle before its next interaction, 
note first that where Cmeias denotes the total inelastic scattering cross-section, 
and N denotes the number of targets per unit volume, the product -/V<7j no i as 
gives the probability per unit path-length of the incoming particle track, that 
the projectile will undergo some type of inelastic interaction. The reciprocal of 
this value therefore gives the mean free path X ine i as for inelastic scattering. In 
other words, Xineias = ^-/N<7 ine i as is the mean distance travelled by a particle 






G 



between inelastic collisions. Similarly, \ e i as = Na e \ as gives the probability per 
unit path length, that the projectile will undergo some type of elastic interaction, 
and the reciprocal 1/Na e i as gives the mean free path for elastic scattering. The 
total mean free path A t is then 

A 



inelas elas 



Now, the path-lengths between interactions will have the following inverse- 
exponential distribution 

P (x) = lexp 

with A t as the mean of the distribution. Given a randomly generated number 
li\ G (0, 1), the distance travelled I by a particle before its next interaction can 
be randomly generated according to 

1 / —x 



The indefinite integral of the path-length distribution is 



F(x) = J -i- exp (^-^—^J dx' = - exp 



— x 

a7 



hence 



IH = J y exp ^— ^ ) (Li 

F(l)-F(0) 



Therefore, 



and 



1 — ji\ = exp 



from which it follows that 



ln(l-/ii) = ^ 

At 



I = -At ln(l - fii) 



Given that \i\ is uniformly distributed in (0,1), 1 — /ii will also be uniformly 
distributed in (0,1), hence, for the purposes of generating random numbers, [i\ 
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can be substituted into this equation in place of 1 — fii, to obtain the final 
expression: 

I = -A t ln(/ii) . 

Given random numbers G (0,1), one uses this expression to generate the 
random path lengths between collisions. 

To decide whether the next collision is to be elastic or inelastic, one can use 
the fact that the probability of an inelastic collision Pi ne ias is specified by 

inelas 

Pineias \— 1 \— 1 

inelas ' elas 

A random number fj,2 is generated, and if /j,2 < Pineias, then the collision is 
chosen to be inelastic, otherwise it is chosen to be elastic. 

The types of possible inelastic collision then depend upon the type of pro- 
jectile under consideration. For example, in the case of an incoming gamma-ray 
photon, there could be a photoelectric interaction, a Compton-effect interac- 
tion, or a pair-production interaction. The relative probabilities of the different 
interaction types depend upon the energy of the incoming particle, and given 
these relative probabilities, a randomly number /13 £ (0, 1) can determine which 
interaction is selected. Let k denote the linear attenuation coefficient of energy 
E gamma rays in the medium, let r denote the photoelectric attenuation coef- 
ficient, and let zu denote the Compton attenuation coefficient. If /i3 is between 
and t/k, then photoelectric absorption can be deemed to occur; if ^3 is be- 
tween t/k and (r + w)/k, then Compton scattering can be deemed to occur; 
and if ^3 is between (r + vo)/k and 1 , then pair production can be deemed to 
occur (Mackie 1990, p552). The energy lost as a result of the interaction is 
chosen, once more, by the generation of random numbers, and by the energy 
loss distributions specific to each type of interaction. 

Given the relevant differential cross-section a(9, </>) for the incoming energy 
E, random numbers can be generated which specify the polar angle and az- 
imuthal angle with which the projectile scatters. Assuming the incoming pro- 
jectile is travelling in the direction # = 0, the probability p(9) of scattering at 
polar angle 9 is 

= cr(6>,</>) sing 
J o Tr cr(6»',0)sin6"rf0' ' 

The integral Jq <j{9' , <fi) sa\9'd9' of the differential cross-section over all polar 
angles gives a quantity with dimensions of Area • Angle -1 , proportional to the 
probability, per unit azimuthal angle, of scattering at any polar angle. The 
expression a{9, <f>) sing gives a quantity, with the same dimensions, proportional 
to the probability of scattering, per unit azimuthal angle, at polar angle 9. The 
ratio of the latter by the former gives the probability p{9) of scattering at polar 
angle 9. Hence, given the selection of a random number [14 G (0, 1), the angle 9 
such that 
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p(6')d6' 



J^a(9',(t>) sm9'd6' 

is selected. The azimuthal angle <j> is determined by a uniform distribution over 
(0, 27r), and the selection of another random number fi 5 e (0, 1). 
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